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Abstract 

A consistent description of neutrino oscillations requires either the quantum- mechanical 
(QM) wave packet approach or a quantum field theoretic (QFT) treatment. We com- 
pare these two approaches to neutrino oscillations and discuss the correspondence 
between them. In particular, we derive expressions for the QM neutrino wave packets 
from QFT and relate the free parameters of the QM framework, in particular the effec- 
tive momentum uncertainty of the neutrino state, to the more fundamental parameters 
of the QFT approach. We include in our discussion the possibilities that some of the 
neutrino's interaction partners are not detected, that the neutrino is produced in the 
decay of an unstable parent particle, and that the overlap of the wave packets of the 
particles involved in the neutrino production (or detection) process is not maximal. 
Finally, we demonstrate how the properly normalized oscillation probabilities can be 
obtained in the QFT framework without an ad hoc normalization procedure employed 
in the QM approach. 
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1 Introduction 



It is well known by now that neutrino oscillations can be consistently described either in 
the quantum-mechanical (QM) wave packet approach, or within a quantum field theoretic 
(QFT) framework.^ In the QM method |T], |5|-P^, neutrinos produced in weak-interaction 
processes are described by propagating wave packets, the spatial length of which is related 
to the momentum uncertainty at neutrino production, and the detected states are also 
described by wave packets, centered at the detection point. The transition (oscillation) am- 
plitude is then obtained by projecting the evolved emitted neutrino state onto the detection 
state. In the QFT treatment |9|-[TT| ^^-|26[, one considers neutrino production, propaga- 



tion and detection as a single process, described by a tree-level Feynman diagram with the 
neutrino in the intermediate state (see fig. |1|). Neutrinos are represented in this framework 
by propagators rather than by wave functions. Both approaches lead to the standard for- 
mula for the probability of neutrino oscillations in vacuum in the case when the decoherence 
effects related to propagation of neutrinos as well as to their production and detection can 
be neglected. They differ, however, in the way they account for possible decoherence effects, 
with the QFT approach leading to a more consistent and accurate description. The QM 
method treats neutrino energy and momentum uncertainties responsible for these effects in a 
simplified way; in addition, it involves an ad hoc normalization procedure for the transition 
amplitude that is not properly justified. 

The goal of the present paper is to compare the two approaches and establish a re- 
lationship between them, as well as to clarify some of the procedures that are employed 
in the QM method from the more general and consistent QFT standpoint. Some work in 
that direction has been done before. In [jl2| neutrino wave packets were derived starting 



from the QFT formalism (see also a discussion in 0). In [^, a comparison of the QM 
and QFT approaches was presented for the special case of Mossbauer neutrinos, i.e. neu- 
trinos produced and detected recoillessly in hypothetical Mossbauer-type experiments (see, 
e.g., p6|, ^ and references therein). The new results obtained in the present work 
include a more advanced and general study of the QFT-based derivation of the neutrino 
wave packets (including the possibility that some of the external particles are not detected), 
matching of the QFT and QM expressions for the neutrino wave packets, study of the gen- 
eral properties of the wave packets describing the neutrino states (including their energy 
uncertainties in the case when neutrinos are produced in decays of unstable particles) and 
clarification of the issue of normalization of the oscillation probabilities in QM and QFT. 

The paper is organized as follows. To make the presentation self-contained, in Sees. ^ 
and 1^ we review, respectively, the QM wave packet formalism and the QFT approach to 
neutrino oscillations. Sections § - P contain our main results. In Sec. § we discuss how the 
neutrino wave packets, which are a necessary ingredient of the QM approach, can be derived 
starting from the QFT formalism. Next, we consider some general properties of the neutrino 



""^ Although in a number of sources a plane wave approach to neutrino oscillations is employed, it is actually 
marred by inconsistencies and, if applied correctly, does not lead to neutrino oscillations at all ^ . 
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wave packets and discuss the conditions under which they can be approximated by Gaussian 
wave packets. Using the case of Gaussian wave packets as an example, we then discuss how 
the QFT-derived wave packets can be represented in the form usually adopted in the QM 
treatment. We also find expressions for the effective parameters describing the QM wave 
packets in terms of the more fundamental input parameters of the QFT framework. Next, we 
discuss the neutrino energy uncertainty in the case when neutrinos are produced in decays of 
unstable particles. In Sec. ^ we consider the problem of normalization of the neutrino wave 
packets in the QM framework and show how the normalization problem is solved in a natural 
way in the QFT-based approach. In Sec. ^ we discuss how one can relax some assumptions 
usually adopted in the QM and QFT approaches. Those include the assumption that the 
maxima of wave packets of all particles involved in the neutrino production (or detection) 
process meet at one space-time point, as well as the assumption that the mean momenta 
of the emitted and detected neutrino wave packets coincide. We summarize our results and 
conclude in Sec. ^ Some technical material is included in Appendices A and B. 



2 Review of the QM wave packet formalism 

We start with some generalities that are common to QFT and QM and then move on to 
review the QM wave packet approach to neutrino oscillations. We shall use the natural 
units h = c = 1 throughout the paper. 

In quantum theory, one-particle states of particles of type A can be written as 

\A) = J[dp]fA{p,P)\A,p), (1) 

where |A, p) is the one-particle momentum eigenstate corresponding to momentum p and 
energy Ea{p) (for free particles, Ea{p) = ^p^ + m^, tua being the mass of the particle), 
/^(p, P) is the momentum distribution function with the mean momentum P, and we use 
the shorthand notation 

[dp] ^ d^P^ (2) 



For particles with spin, the states \A) and |v4,p) depend also on a spin variable, which we 
suppress to simplify the notation. We will also often omit the second argument of Ja where 
this cannot cause confusion. 

We choose the Lorentz invariant normalization condition for the plane wave states |v4, p): 

{A, p'\A, p) = 2EAip) (27r)35(3)(p - p') . (3) 
The standard normalization of the states = 1 then implies 

IfslWp)!^-!. W 
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The quantity ■\/2£'a(p)/a(p) is actually the momentum-space wave function of A: ^2i?^(p) ) 
/^(p) = {p\A). The time dependent wave function is ^2EAip)fAip)e~'^^^^^^ = {p\A{t)), 
where \A(t)) = e~^^^\A) and H is the free Hamiltonian of A. The coordinate-space wave 
function \l/^(t,x) is the Lorentz-invariant Fourier transform of (p|y4(t)):^ 



or 



vl/^(t, x) = 1 H /^(p)e-^^-(P)*+^P^ . (6) 

In the QFT framework, it can be written as 

^a{x) = {0\^a{x)\A), (7) 

where x = {t, x) and "^aIx) is the second-quantized field operator of A. Using the standard 
decomposition of the field "^aIx) in terms of production and annihilation operators, one 
can readily obtain (|^) from (|^ and (|1]). Note that expressions (^) and (|^) can describe 
both bound states and propagating wave packets (in the case of bound states or particles 
propagating in a potential, the relation -^'^(p) = a/p^ + fTi\ simply has to be replaced by 
the proper dispersion relation). A wave packet is obtained when the momentum distribution 
function /a(p, P) is sharply peaked at or close to a nonzero mean momentum P,'^ i.e. when 
the momentum dispersion ap satisfies ap <^ |P|; for the rest of this section we will assume 
this to be the case. The wave function (P) then describes a wave packet whose maximum 
of amplitude is located at x = at time t = 0. A wave packet that is peaked at coordinate 
xq at time to is obtained by acting on the state \A) by the space-time translation operator 
giPxo^ where P'^ is the 4-momentum operator. For the coordinate-space wave function this 
yields 

^^(x) = J [dp] /^(p)e-^^^(p)(*-*o)+^P(''-''o) (8) 

instead of eq. (^. 

Eqs. (^) and represent wave packets that propagate with the group velocity v = 
^^dp^^ I p=p = -E^p) ^^"i general spread with time both in the longitudinal direction and 
in the directions transverse to their mean momentum. The spreading is due to the fact 
that different momentum components of the wave packet have slightly different velocities 
p/^a(p). 



^Recall that the Fourier transformation is based on the completeness condition for 1-particlc momentum 
eigenstates, which for our normahzation convention reads / (27r)4g4(p) IP)(Pl ~ ^ipart. (see, e.g., Q, 
eq. (2.39)). Here the right hand side is the unit operator in the subspacc of 1-particlc states and zero in the 
rest of the Hilbert space. Note that the integration measure d^p/ Ea{p) is Lorcntz invariant. 

"^The peak momentum coincides with the mean momentum for symmetric wave packets. 
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Let us now consider neutrino oscillations in the framework of the QM wave packet 
formalism, sometimes also called the "intermediate wave packet" approach. Neutrinos pro- 
duced or absorbed in charged- current weak interaction processes are considered to be flavour 
eigenstates Ua {a = e,fi,T), which are coherent linear superpositions of mass eigenstates uj 
(j = 1, 2, 3) with coefficients given by the elements of the leptonic mixing matrix Uaj- The 
mass eigenstates are represented by the corresponding wave packets. If a neutrino of flavour 
a was produced at time tp at a source centered at xp, its momentum-space wave function 
at a time t > tp is 

{pWMt)) = Yl U:,V'2EAp)fMP^ P) e-^^(P)(*-*-)-P-- . (9) 
j 

Here the subscript P shows that the wave packet corresponds to a neutrino produced at the 
source. Note that the index a at t'ap(t) simply indicates that the emitted neutrino was of 
flavour a at its production time t = tp; it is, of course, no longer so for t > tp. The shape of 
the wave packet of the jth mass-eigenstate neutrino is given by the momentum distribution 
function fjp, which is determined by the mechanism and conditions of neutrino production. 
In the QM framework, however, the neutrino production and detection processes are not 
explicitly taken into account; therefore the functions fjp are postulated rather than deter- 
mined, with the corresponding momentum widths estimated from the localization properties 
of the production process. Usually, the wave packets are taken to be of the Gaussian form 

(^pp J L "^^pp ^ 

where app characterizes the momentum uncertainty of the produced neutrino state, and 
similarly for the state of the detected neutrino. The advantage of Gaussian wave packets is 
that they allow most calculations to be done analytically (the same is also true for Lorentzian 
wave packets, see ref. [^). 

The state of the detected neutrino is described by a wave packet peaked at the 
detection coordinate xp). In the momentum-space representation it is given by 

{pW^d) = J2u;,,/2EAp) fkD{p,P')e-'^-- , (11) 

k 

where the subscript D stands for detection. The momentum distribution functions f^D are 
governed by the properties of the detection process; however, just as for neutrino production, 
in the QM approach these functions are postulated rather than determined. Note that, 
although the assumption P = P' is adopted in most studies, in general there is no reason 
to expect the mean momenta of the produced and detected wave packets to coincide. We 
will discuss this point in more detail in Sec. |[ 

The amplitude for the transition Ua — I'p is obtained by projecting the evolved neutrino 
production state onto the detection state: 

A/3(T,L) = (z/^p|z/„p(tp)), (12) 
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where to is the detection time, T = to — tp > and L = x/? — xp. Performing the projection 
in momentum space, we obtain from (j^) and (|TT]) ^ 

= E ^^.^/^. / 0j-s /.Hp, P) /;d(p, P') e-^^(P)^+-- . (13) 

For future reference, we shall also write this as a superposition of the amplitudes corre- 
sponding to the contributions of different neutrino mass eigenstates: 

Aa^iT, L) = J2 K,Up,MT, L) (14) 
i 

with 

MT. L) = / fAv. P) fUv. P') e-^^^(P)^+^P^ . (15) 

The oscillation probability is given by the squared modulus of the transition amplitude: 
Pq,^(T, L) = P{ua — )• iyfs,T,'L) = |^q,^(T, L)p. Since in most experiments the neutrino 
emission and detection times are not measured, the standard procedure is then to integrate 
P{^a ^13, T, L) over T. This gives 

P{ua ^ up^L) = Paf,{-L) = J dT\A^p{T,L)\\ (16) 

Substituting here the transition amplitude (0) yields, up to a normalization factor, the 
standard probability of neutrino oscillations in vacuum provided that all decoherence effects 
are negligible. The normalization factor can then be fixed by requiring that the oscillation 
probability satisfy the unitarity condition (see Sec. || for a more detailed discussion). 



3 Neutrino oscillations in QFT 

In the QFT approach (which is sometimes also called the "external wave packet" formal- 
ism), neutrino production, propagation, and detection are considered as a single process, 
described by the Feynman diagram shown in fig. |I], with the neutrino in the intermediate 
state. In our overview of the QFT formalism we will mostly follow ref. 0. Assume that the 
neutrino production process involves one initial state and one final state particle (besides 
the neutrino). Likewise, we will assume that the detection process involves only one particle 
besides the neutrino in the initial state and one particle in the final state. The general- 
ization to the case of an arbitrary number of particles involved in the neutrino production 

^Projection in momentum space will turn out to be convenient for our subsequent discussion. The 
coordinate-space projection ^q^(T,L) = / cPx {vaol^) {Mi^ap{tD)) yields, of course, the same result. 
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Figure 1: Feynman diagram describing neutrino production, propagation and detection as 
a single process. 

and detection processes is straightforward and would just complicate the formulas without 
providing further physical insight. All external particles will be assumed to be on their 
respective mass shells.^ 

The states describing the particles accompanying neutrino production and detection 
("external particles") can be represented in the form (|T|). For the initial and final states at 
neutrino production we can write 

\P^) = J [dq] /«(q, Q) li'., q) , \Pf) = j [dk] /p/(k, K) \Pj, k) , (17) 
and similarly for the states accompanying neutrino detection: 

I A) = j m /D.(q , QO I A, q ) , 1^/) = / [^^1 /^/(k'. K') \Df, k') . (18) 

We assume these states to fulfill the normalization condition (^. Some (or all) of the mean 
momenta of the external particles Q, K, Q' and K' may vanish, i.e. the states in eqs. ( [T7| ) 
and (|18D can describe bound states at rest as well as wave packets. 

The amplitude of the neutrino production - propagation - detection process is given by 
the matrix element 

iAai5 = {PfDf\fe^^ - i j d^xUiix) -l|PiA), (19) 

where T is the time ordering operator and 'Hi{x) is the charged- current weak interaction 
Hamiltonian.^ Note that no neutrino flavour eigenstates have to be introduced in the QFT 



^Since only one particle is assumed to be in the initial state of the production process, it must be 
unstable. This will be of no importance for us here because, as was already mentioned, the results are 
easily generalized to the case of an arbitrary number of external particles. Possible instability of the parent 



particle will be taken into account in Sec. 4.4. 

®We consider neutrino production and detection at energies well below the M^-boson mass, so that T-Li is 
the effective 4-fermion Hamiltonian of weak interactions. 
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framework, and the indices a and /5 simply refer here to the flavour of the charged leptons 
participating in the production and detection processes. 

From eq. (19) it is easy to calculate the transition amplitude in the lowest nontrivial (i.e. 



second) order in "H/ using the standard QFT methods. The resulting expression corresponds 
to the Feynman diagram of fig. |l| and can be written as 

lAo^f, = ^2^:^^^, J[dq]fp,{q,Q) j[dk]rpf{k,K) 

X I [dq'] fmiq!, QO I [dk'] fhfik', K') z^^---(g, k; q' , k') . (20) 

Here the sum runs over all intermediate states (i.e. different neutrino mass eigenstates), 
and the quantity A'j'^'^q, k; q', k') is the plane- wave amplitude of the process with the jth 
neutrino mass eigenstate propagating between the source and the detector: 

tAf'"-{q,k;q',k') = j d^xij rf^xs MD(g , A;') e'^^^'-'^'^^"^-"^) 

X I I J^"^' e-ip{-2-x,) . ^ / ^^ g-i(g-fc)(xi-xp) ^ .21) 

J {2ttY - m^- + le ' ^ ' 

Here x\ and x^ are the 4-coordinates of the neutrino production and detection points, the 
quantities Mp{q, k) and Mniq', k') are the plane-wave amplitudes of the processes Pi ^ Pf + 
Uj and Di + Vj — >■ -D/, respectively, with the neutrino spinors Uj{j), s) and Mj(p, s) excluded. 
The choice of the 4-coordinate dependent phase factors corresponds to the assumption that 
the peaks of the wave packets of particles involved in the production process are all located at 
xi = xp at the time ti = tp, whereas for the detection process the corresponding peaks are 
all situated at X2 = x^^ at the time t2 = tp) (we will discuss in Sec. ^how this assumption can 
be relaxed). The integral in the second line of eq. (|2T|) gives the coordinate-space propagator 
of the jth neutrino mass eigenstate. 

It is convenient to switch to shifted 4-coordinate variables x[, x'^^ defined according to 
X\ = xp + x[, X2 = Xp, +X2- Taking into account that 'f + rrij = Uj{p, s)uj{p, s), one can 
then rewrite eq. (|2T|) as 



/d'^p Q-ip{^D-xp) r 
TTT-TT^ I yWoMjp{q, k) e 



X y rfVa V2poM,^(g,A;')e"*^'^+P""^"2, (22) 
where 

M,p(g, k) ^ ^^^Mp{q, k) and M,^(g', k') ^ Mp{q', k')"^^ (23) 
V ^Po V ^Po 
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are the full amplitudes (with the neutrino spinors included) of the processes Pj — )■ P/ + Vj 
and Di + Vj Df, respectively, and we have taken into account that the matrix elements 
Mjp{q, k) and MjD{q', k') involve the left-handed chirality projection, so that only the left- 
handed spinors Ujiip) and Ujiip) contribute to the sum over the neutrino spin variable s. 

Substituting (|2^ into eq. (|20|) , we finally obtain 

^A^, = ^ E J ^ $,p(p°, p)$,M/, p) . (24) 

Here the so-called overlap functions $jp(p°,p) and ^jr,{p^,v) ^-^e defined as 

$ip(/, V) =j d'x[e'^< I [dq] j [dk] /p.(q, Q) /^/k, K) e"^(^-'=)<M,p(g, k) , (25) 

$,d(/,p) = j d^e-'^^'- j[dq'] y"[rffc']/z,.(q',QO/^,/(k',K')e-^(^'-'=>^M,z,(g',A;'). 

Note that they are independent of xp and xd. Expressions (^) and ( ^5] ) are the main 
results of the QFT-based approach to neutrino oscillations |0, |12| . 

4 Comparing the QM and QFT approaches to neu- 
trino oscillations 

Let us now compare the results of the QM and QFT approaches to neutrino oscillations. 
Consider first the transition amplitude (p^ ) obtained in the QFT formalism. The integration 
over the neutrino 4-momentum in this expression can be done in different order. Here it 
will be more convenient for us to integrate first over p^ and then over p (the opposite order 
will be used in Sec. |^). Since the distance L between the neutrino source and detector is 
macroscopic, the phase factor in the integrand of eq. (|24D undergoes fast oscillations and 
the integral is strongly suppressed except when the intermediate neutrino is on the mass 
shell. Thus, the dominant contribution to the integral is given by the residue at the pole of 
the neutrino propagator at p^ = Ej{p) — iej where 

E,ip) = ^p^ + m]. (26) 
Eq. (^) can therefore be rewritten as 

= e(T) J2 U:,U,, I ^ $,p(i?,(p), P)^ME,{P), p) e-^^^(P)^+^P^ . (27) 

where 6(x) is the Heaviside step function. 

''The contribution of the pole at = —Ej{p) + ie is strongly suppressed due to an approximate conser- 
vation of mean energies at production and the fact that i?pi(Q) — £'p/(K) > 0. 



9 



4.1 Deriving neutrino wave packets in the QFT-based approach 



Let us now compare eqs. (^) and (|TB|). We see that the two equations are of the same form 
and actually coincide if we identify the QM wave packets as 

fjpip) = <l>,p(E,(p),p) , fjDip) = <^*DiEjip),p) , (28) 

where the functions and were defined in eq. (p5|). 

The obtained result can be easily understood. Indeed, as follows from the definition 
of $jp(p°,p), for p° = Ej{p) (i.e. on the mass shell of Uj) this quantity is the probability 
amplitude of the production process in which the jth mass eigenstate neutrino is emitted 
with momentum p; but this is nothing but the momentum distribution function of the 
produced neutrino, i.e. the momentum-state wave packet fjp{p). A similar argument applies 
to the neutrino detection process and fjD^p). The wave packets /jp(p) and fjoip) in 
eq. ( |28D are not normalized according to (H), though they can be easily modified to satisfy 
this condition. However, as we shall see in Sec. ^, this is not necessary and actually would 
be misleading. 

An alternative method of deriving neutrino wave packets in the QFT framework, based 



on the S-matrix approach, was suggested in [|12|; the obtained results are equivalent to those 
in eqs. and (^). 

Let us now consider the wave packet describing the produced neutrino state in more 
detail (the state of the detected neutrino can be studied quite analogously). According to 
(|28|) , the momentum distribution function fjp{p) characterizing the state of the emitted 
neutrino of mass rrij is essentially given by the on-shell function $jp(£'j(p), p). Since the 
matrix element Mjp{q, k) is a smooth function of the on-shell 4-momenta p and q, whereas 
the wave packets of the external states are assumed to be sharply peaked at or near the 
corresponding mean momenta, one can replace Mjp by its value at the mean momenta and 
pull it out of the integral. Eqs. (p5[) and (|28|) then yield 



/,p(p) ^ M,p(g, K) J d'x e^^.(p)*-P-y [dq] J [dk] /p,(q, Q) /^/k, K) e-'^^-'^^ , (29) 

where the 4-momenta Q and K are defined as 

g = (Ep,(Q), Q), K = {Epf{K),K). (30) 

From eq. (|29| ) (or eqs. ( pSj) and (pH])) one can draw some important conclusions about the 
properties of the neutrino momentum distribution functions fjp{p) which determine the 
emitted neutrino wave packets: 

• Since the quantities /pi(q, Q) and /p/(k, K) depend only on the properties of the ex- 
ternal particles, and the j-dependence of the matrix elements Mjp{p, q) comes through 
the on-shell neutrino spinor factors [(2p°)~^/^Uj(p, s)]pO^Ej{p)^ which depend on j only 
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through the neutrino energy, the functions /jp(p) depend on the index j solely through 
the neutrino energy Ej{p). This, in particular, means that for ultra-relativistic or 
quasi-degenerate in mass neutrinos the momentum distribution functions of all neu- 
trino mass eigenstates are practically the same (provided that their energy differences 
\Ej — Ek\ are small compared to the energy uncertainty (Xep). 



Because the integral over the 3-coordinate x in eq. (p9| ) yields S^^^ (q — k — p), and the 
momentum distribution functions /pi (q,Q) and /p/(k, K) are sharply peaked at or 
near their respective mean momenta Q and K, the neutrino momentum distribution 
functions /,p(p) are sharply peaked at or close to the momentum P = Q — K, with the 
width of the peak app dominated by the largest between the momentum uncertainties 
of the states of Pi and Pf. 



Taking into account eq. (^, eq. 



can be rewritten as 



pO=E,(p) 



(31) 



Thus, the momentum distribution function that determines the wave packet of the emitted 
neutrino is essentially the 4-dimensional Fourier transform of the product of the coordinate- 
space wave functions of the external particles participating in the neutrino production pro- 
cess, taken under the condition that the components of the neutrino 4-momentum are on 
the mass shell. Eq. (|3TD can be readily generalized to the case when more than two external 
particles participate in the neutrino production process: the expression \Efp.(t, x)\E'p^(t, x) 
in the integrand of (|3ll) should simply be replaced by the product of the wave functions of 
all particles in the initial state of the production process and of complex conjugates of the 
wave functions of all particles in the final state (except the neutrino). 

The neutrino wave packet in coordinate space il)jp{x) is obtained from eq. (|31| ) by 
performing the Fourier transformation over the 3-momentum variable p according to the 
transformation law (|^), which gives 

i,^p{x)c^u,{P,s)M,p{Q,K) ! d^'x'^pXx'Wpix') 



e(t - 1') 



|x-x'| 2(27r)2 7o I 2 + ^2 



g-iA/p2+m2(t-t') |-^ip|x-x'| _ g-ip|x-x'|j _ |-22^ 



Here p = |p| and we have used eq. (|23|) to extract the p-dependent factor (2p°)~^/^ = 
2~^/^(p^ + nij)~~^^^ from Mjp{Q, K). The integral over p in eq. (^) can be expressed in 
terms of the modified Bessel function Ki [0, giving 



(2x) 



\/—{x — x'Y 



(33) 
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where (x — x')^ = {t — t'Y — (x — x')^. The integral in the second hne of eq. (|3^) is greatly 
simplified in the limit of vanishing neutrino mass: 



Note, however, that in this limit all neutrino species travel with the same speed and therefore 
the wave functions in eq. ( p4D cannot describe decoherence due to the separation of wave 
packets. In order to take possible wave packet separation effects into account the more 
accurate expression ( P^ has to be used. Alternatively, one can employ the momentum- 
representation wave function (|31|). 



Expression (|33D for the wave function of the produced neutrino state \&jp(t, x) allows 



a simple interpretation. Note that mjQ{t — t')Ki{mj\/ —{x — x'Y)/ \/—{x — x' )^ is 
the scalar retarded propagator in the coordinate representation. Therefore the neutrino 
wave packet (^) is essentially the convolution of the neutrino source (the role of which is 
played by the neutrino production amplitude \Efp^(x)M,p\l/p^(x)) with the retarded neutrino 
propagator, in full agreement with the well known result of QFT. Note that only the scalar 
part of the propagator contributes to \E'jp(x); this is because the coordinate space and 
momentum space neutrino wave functions \l/jp(a;) and /,p(p) are scalars in our formalism. 
The spinor factors are included in the matrix elements Mjp and Mju (note that these 
quantities are also scalar, whereas the amputated matrix elements Mjp and Mjo, i.e. those 
with the neutrino spinors removed, have spinorial indices). 

In our discussion of the wave packets of the emitted neutrino states, we were assuming 
that the momentum distribution functions of all the external particles accompanying neu- 
trino production are known. This implies, in particular, that all particles in the final state 
of the production process are "measured" , either by direct detection or through their inter- 
action with the medium in the process of neutrino production. It is quite possible, however, 
that some of the particles accompanying neutrino production escape undetected; this is, e.g., 
the case for atmospheric or accelerator neutrinos born in the process vr^ — )■ /x^ -|- i^nii'fj,), 
in which the final state muon is normally not detected. It is also possible that some of the 
particles accompanying neutrino detection are "unmeasured". How can one determine the 
neutrino wave packets in those cases? 

To answer this question, let us recall that the momentum uncertainty o"pp characterizing 
the emitted neutrino depends in general on the momentum uncertainties of all the external 
particles at neutrino production and is dominated by the largest among them (see the 
discussion after eqs. (|^) and (|30|) ). In particular, in the case of Gaussian wave packets, one 
has P, ra 



(^lp = ^lpi + ^lpf (35) 

For more than two external particles at production, the sum on the right-hand side of 
this relation would contain the contributions of the squared momentum uncertainties of all 
these particles. Now, if a particle goes "unmeasured" in the neutrino production process. 
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its momentum uncertainty cannot affect the momentum uncertainty of the emitted neu- 
trino state and therefore can be neglected. To put it differently, undetected particles are 
completely delocalized, and therefore, according to Heisenberg's uncertainty relation, have 
vanishing momentum uncertainty. This means that undetected particles can be represented 
by states of definite momenta, i.e. by plane waves. If, for example, the particle Pf at pro- 
duction is undetected, one has to replace in eq. (pOf ) the momentum distribution function 
/pj(k, K) by [(27r)^/\A^]5''^^(k — K) where V is the normalization volume, and in eq. (|3T|) 
the coordinate-space wave function \l/pj(x) by e~*^^/y^^!Ep/(K)V, with eqs. (^) - (|3^ ) 
modified accordingly. The mean momentum of the neutrino state depends, of course, on 
the momentum of the undetected particle; if the latter can take values in some range, the 
same will be true for the mean momentum of the emitted neutrino state. In this case the 
flux of emitted neutrinos will be characterized by a continuous spectrum. 

In most of our discussion in this subsection we concentrated on the wave packets of the 
produced neutrino states. Our consideration, however, applies practically without changes 
to the detected neutrino states; the corresponding formulas can be obtained from eqs. 
and (|3TD -(|34D by replacing, where appropriate, e^^^ — >■ 6"*^"^, \l/pj — "^^i and \E'p/ — > "^dj- 

4.2 General properties of neutrino wave packets 

We have already considered some of the general properties of the neutrino wave packets in the 
previous subsection. In particular, we have found that the momentum distribution functions 
fjp{p) of mass-eigenstate neutrinos uj depend on the index j only through the neutrino 
energy Ej{p), and that the functions fjp{p) are sharply peaked at or near the momentum 
P = Q — K, with the width of the peak dominated by the largest between the widths 
of the functions fpi and fpf. Further insight into the general properties of the neutrino 
wave packets can be gained by comparing expressions dl) with their pi ane-wave limits. 
If the external particles were described by plane waves, the quantities ^jp{Ej{p),p) and 
^j£){Ej{p), p), which determine the neutrino wave packets, would be just equal to the matrix 
elements of the neutrino production or detection processes divided by the factor \/2EV for 
each external particle and multiplied, correspondingly, by {27r)^S^'^\q — k^p). The latter 
factors represent energy-momentum conservation at the production and detection vertices. 
As follows from (^5]), in the case when the external particles are described by wave packets, 
the quantities $jp(_Ej(p), p) and ^jD{Ej{p),p) (and therefore the momentum distribution 
functions of the neutrino wave packets) correspond to "smeared (5-functions" , representing 
approximate conservation of the mean energies and mean momenta of the participating 
particles. How exactly this smearing occurs will depend on the form of the wave packets of 
the external particles, and to move ahead one has to specify this form. 

A particularly useful and illuminating example of a specific form of the external wave 
packets, and the one most often used in the literature, is the case of Gaussian wave packets. 
We will employ this example to illustrate the general properties of the neutrino wave packets. 
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Let us discuss first the conditions under which an arbitrary wave packet can be accu- 
rately approximated by a Gaussian one. For simphcity, we will consider here 1-dimensional 
wave packets. This is a good approximation in the case when the distance between the 
neutrino source and detector is very large compared to their sizes, so that the neutrino mo- 
mentum is practically coUinear with L = x/) — xp (the generalization to the 3-dimensional 
case is straightforward). Consider a wave packet described by a momentum distribution 
function f{p), sharply peaked at some value Pq of the momentum. We can write this 
function in the exponential form as 

f{p) = e-f(p) , where g{p) = -\n[f{p)] . (36) 

The Gaussian approximation corresponds to the case when the integral over p of the function 
f{p) multiplied by any function of p that is smooth in the vicinity of Pq can be evaluated 
in the saddle point approximation. Indeed, in this approach one expands the function g{p) 
around its minimum at p = Pq a-nd keeps the terms up to and including the quadratic one: 

gip) ^ giPo) + i g"iPo)ip - Pof • (37) 

This precisely means the wave packet f{p) is approximated by the Gaussian one. The 
validity condition for this approximation is given in terms of the derivatives of the function 
g at Pq: 

1 , 1 



\9 



(^^)(Po)|«7T|/(Po)r. (3^ 



4! ' 2 

It can be satisfied for a wide range of functions f{p)- However, it is easy to construct wave 
packets for which it is not satisfied. Consider, e.g. a class of wave packets 

with integer n and C„ a constant, which can be found from the normalization condition for 
/(p). It is easy to check that condition (^) is equivalent to l/4n <C 1. Thus, the momentum 
distribution functions (^) can be accurately approximated by the Gaussian ones only when 
n ^ 1. This condition, in particular, is not satisfied for Lorentzian wave packets. 



4.3 Matching the QFT and QM neutrino wave packets 

Let us now discuss how one can match the QFT and QM wave packets of neutrinos. Using 
the case of Gaussian wave packets as an example, we shall find out how the effective param- 
eters describing the QM wave packets can be expressed in terms of the more fundamental 
parameters entering into the QFT approach. 

We start by introducing some notation (we mostly follow 0, [1^ here). The coordinate 
uncertainty a^pi characterizing the wave function of the initial state particle Pi in the 
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neutrino production process is related to its momentum uncertainty apPi by 



o'xPiO-pPi = - , (40) 

and similarly for all other external particles. One can also introduce the effective coordinate 
uncertainty of the production process a^p, which is connected to the effective momentum 
uncertainty of this process Cpp defined in eq. (^) by a relation similar to (^Op, or equivalently 

1 1 1 

T^ = T^ + T^- (41) 

'^xP ^xPi ^xPf 

This formula has a simple physical interpretation: since the neutrino production process 
requires an overlap of the wave functions of all the participating particles, the effective 
uncertainty of the coordinate of the production point is determined by the particle with 
the smallest coordinate uncertainty. This is in accord with the already discussed fact that 
the effective momentum uncertainty at production o"pp, which determines the momentum 
uncertainty of the produced neutrino, is dominated by the largest among the momentum 
uncertainties of all the external particles involved in neutrino production. 

Next, we define the effective velocity of the neutrino source vp and its effective squared 
velocity Sp as 

- ^2 ( ^Pi ,^Pf\ ^ -2 ( V|>, ^Pf \ , . 

\'-'xPi ^xPf J \^xPi ^xPf J 

If vpi ~ vpj, they are approximately equal to, respectively, the velocity and squared velocity 
of the particle with the smallest coordinate uncertainty. We will also need the quantity (XeP 
defined through 

^Ip = ^Ip{^p - Vp) = cr^p Xp . (43) 

This quantity can be interpreted as the effective energy uncertainty at neutrino produc- 
tion 0. It can be also shown that < Ap < 1, i.e. < a^p < Cpp. 

We can now discuss the results obtained in the QFT framework in the case when the 
external particles are represented by Gaussian wave packets. The function <l>jp(£'j(p), p), 
which coincides with the momentum distribution function of the emitted mass-eigenstate 



neutrino Vj, can be written as 0, |T2 



$,p(E,(p), ^)=NpMjp{Q,K) ^ exp [-(?p(E,(p), p)] , (44) 

ePO^pP 



where 



(2vra^pJ3/4(2,ra^p )3/4[2i^^^(Q).2Ep,(K)]i/2 ^^^^ 
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is the normahzation factor and 



9p[Ej[p), p) - —-2 + T-2 ■ l4bj 

Here 

P = Q-K, Ep = Ep,iQ)-Epf{K), (47) 

and Ej{p) was defined in eq. (pB]). Note that in the hmit when the external particles are 
represented by plane waves (cTppj — )■ 0, appj 0), the first equation in ( P^D yields 

<l>,p(^,(p), p) = (27r)^5[Ep,(Q) - Epf{K) - E,{P)] 6^^\Q - K - P) 

(48) 

^2Ep,{Q)V-2Epf{K)V 



as discussed in the previous subsection. From eq. ( p5D and the fact that (XeP < cTpP it 
follows that in this limit the momentum uncertainty app and energy uncertainty aep of the 
produced neutrino state vanish as well; thus, if the external particles are described by plane 
waves, then so is the produced neutrino. As follows from eq. (0), the plane wave limit 
corresponds to exact energy and momentum conservation at production.® This can also be 
seen from eqs. (|4^) and (|46|): indeed, in the limit cxep, Upp — )■ the right hand side of (Q) is 
proportional to the product of the energy and momentum conserving (5-functions. For finite 
values of app and a^p, eqs. (^) and (^61) yield Gaussian-type "smeared delta functions", 
i.e. describe approximate conservation laws for the mean momenta and mean energies of 
the wave packets, for which are responsible, respectively, the first term and the second term 
in (M). 



Let us now try to cast expressions (|^) and ( ^q) into the form usually adopted in the 
QM wave packet approach to neutrino oscillations. We want to reduce $jp(i?j(p), p) to an 
expression similar to that in eq. (|10|). To this end, we expand the neutrino energy around 
the point p = P and keep terms up to the second order: 

E,{p) ^ E, + v,(p - P) + ^(5^'' - v^v^^)ip - P)\p - Py . (49) 



Here 



E,.E,(P), v,.^ 



P=p Ej 



dp 

The lower index j corresponds to the neutrino mass eigenstates, while the upper indices k 
and / number the components of the 3- vectors p, P and Vj (i.e. Vj is the kth component of 



^ Since the neutrino energy and momentum are completely determined by those of the external particles in 
this case, only one neutrino mass eigenstate can be produced in any given interaction process, and therefore 
no oscillations are possible in the plane wave limit H . 
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the group velocity of the j'th neutrino mass eigenstate). The function gp{Ej{p),p) defined 
in ( ^BD can then be written as 

gpiE.ip), p) = (p - P)' a'' {p - Py - (5\p - Pf + 7, , (51) 



where 



Xp 5'^ + (v, - vpf (v, - vpY + ^L_^^s'i _ 



(52) 



2 



P'' = -—^iE,-Ep){v,-vpr, 7, = ^ \. ■ (53) 
Let US now try to represent gp{Ej{p), p) in the form 

gp{Ej{p), p) = {p- Peff)' (P - Pcff)' + 7i (54) 

with 

Pcfr = P + <5, (55) 

where the parameters 6 and are to be determined by comparing eqs. (0) and (|5lD . 
They describe, respectively, a shift of the neutrino mean momentum compared to the naive 
expectation (p) = P and a modification of the overall normalization of the neutrino wave 
function. The effective momentum uncertainty characterizing the QM neutrino wave packet 
can be obtained by diagonalizing the matrix a. The squared uncertainties of the different 
components of the neutrino momentum are given, up to the factor 1/4, by the reciprocals of 
the eigenvalues of the matrix a. In general, these eigenvalues are different, i.e. the neutrino 
momentum uncertainty is anisotropic. This actually means that expression (|TUp for a 3- 
dimensional Gaussian wave packet is oversimplified: its exponent has to be replaced by 

-[(p^ - p^sfK^^P? + ip' - pi^?/^{^ipf + {'P' - p:k?/^{^;p?]- 

Comparing eqs. ( ^Ij) and (|5T|) , one finds that the shift 5 of the neutrino mean momentum 
satisfies the equation 

2 a^^ 5' = /3^ (56) 

whereas the parameter 7^ is given by 

7,=7,-5'=a''=5' = 7,-^5'=/3^ (57) 

The full solutions of eqs. (|56[) and ( ^Tl) are given in Appendix A; here we present the results 
obtained in the leading order in the small parameter {Ej — Ep)/Ej.^ The diagonalization 
of the matrix a*^' gives in this limit 

(^x ^2 _ ( y ^2 _ 2 ^ _ J_ (Vj " Vp)^ . . 

\"pPeS) ^pP "eP 



^That this parameter is indeed small can be readily seen from eqs. (^) and (^6|): the function is 
strongly suppressed unless \Ej — Ep\ < max{(Tep, cfpp} <C Ej. 
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where the z axis was chosen in the direction of Vj — vp. For b and 7j one finds 

{E, - Ep){v, - vp)' ^ _ {E, - Epf Xp 



Xp + (vj - vp)2 ' ^ Aalp Xp + (vj - Vp; 

From eq. (|58D it follows that the neutrino momentum uncertainty in the direction of v^ — vp is 
smaller than those in the orthogonal directions. Note that for non-relativistic sources (up -C 
1) the direction of — vp essentially coincides with that of the mean neutrino momentum, 
and eq. ( ^8]) means that the longitudinal uncertainty of the neutrino momentum is smaller 
than the transversal ones. To understand this property qualitatively, one can imagine the 
neutrino production region (i.e. the region where the wave packets of the particles involved 
in the production process have significant overlap) to be approximately spherical with radius 
of order a^p. Then, the transverse extent of the neutrino wave packet will also be 0{axp)- 
On the other hand, its longitudinal spread is determined by the duration of the production 
process (i.e. the time interval during which the wave packets have significant overlap), which 
is given by ~ a^p/Sv, where 6v is the relative velocity of the two external particles. 

Summarizing the results of the current subsection, we can see that the QM neutrino 
wave packets can match those obtained in the QFT framework if one applies the following 
changes to the QM results: 

• The momentum uncertainties of the neutrino mass eigenstates are replaced by the 
effective ones, defined in eq. (|58|) . They are in general different in different directions. 

• The mean momentum P is shifted according to P — Peg = P + <5, where the compo- 



nents of 5 are given in eq. (^9]). 



• The wave packet of each neutrino mass eigenstate gets an extra factor Nj = exp[— 7^], 
where 7^- is given in (^9]). 

From the last point one can see that if the differences of the energies of different neutrino 
mass eigenstates are small compared to the energy uncertainty aep, 

Ayr?.?- 

I^i-^.l^^ « ^eP, (60) 

the additional factors Nj are essentially the same for the wave functions of all neutrino mass 
eigenstates and can be included in their common normalization factor. If, on the contrary, 
this condition is violated, the coherence of the emission of different neutrino mass eigenstates 



will be lost |T5|. 



As follows from eq. (|55D, the effective momentum uncertainties that should be used to 
describe the wave packets of emitted neutrinos in the QM formalism are not just equal to 
the true momentum uncertainty at production cr^p, as naively expected, but also depend on 
the energy uncertainty a^p, which is an independent parameter, as well as on the neutrino 



18 



velocity Vj and the effective velocity of the neutrino production region vp. Except for Vj ~ 
vp, the momentum uncertainty along the direction of Vj — vp is dominated by the smaller 
between app and Cep, which turns out to be cxep. This is related to the fact that neutrinos 
propagate over macroscopic distances and therefore are on their mass shell, which enforces 
the relation Ej{p)aep — |p|o"pPefr (see Sec. 5.2 of ref. 



In the limit cXep, vp — j- 0, which 
corresponds to a stationary neutrino source approximation the effective longitudinal 
momentum uncertainty u^p^^ vanishes, even though the true momentum uncertainty cipp 
is nonzero. This implies an infinite coherence length, in accordance with the well known 
result for the stationary case B EO]. It also confirms the expectation that wave packets 



of Mossbauer neutrinos, which are emitted in a quasi-stationary process, have a very large 
spatial extent . 

We have discussed here the matching of the QFT and QM wave packets of the produced 
neutrino states; for the wave packets of the detected states the consideration is completely 
analogous. 



4.4 The case of an unstable neutrino source 

Let us now consider the situation when neutrinos are produced in decays of unstable parti- 
cles. Once again, we will assume that the external particles are described by Gaussian wave 
packets. Compared to the standard formalism that leads to eqs. ( ^l)) and (|^), one now has 
to introduce the following modifications: 

1. The energy of the parent particle Pi acquires an imaginary part, i.e. one has to replace 
Epi{ai} ^P»(q) - zr/2, where T = [mp,/ Ep,{ci)]To ~ [mpi/ Epi{Cl)]To, Tq being 
the rest-frame decay width of Pi. This amounts to replacing the energy difference Ep 
defined in (^) according to Ep — )• Ep — iT/2. 

2. The integration over time in the formula for $jp in eq. (^) now has to be performed 
from to oo rather than from — oo to oo (assuming that t = is the production time 
of the parent particle Pi). 

As a result of these modifications, eqs. (0), (0) get replaced by 

$,p(E,(p), p) = Np M,p{Q, K) -jl-^ exp [ - ^^^P^] h , (61) 



where 



with 



/•oo nOG 

Jo J~tp 

a = alp, 6 = Ej(p) -Ep-vp(p-P), b = b + iV/2, (63) 
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and tp > being the time of maximum wave packet overlap in the production region (see 
Sec. ^). The p-dependent Gaussian factor in eq. (0) describes, as before, an approximate 
conservation of mean momenta at neutrino production, whereas the factor Ii is responsible 
for an approximate conservation of mean energies. Unlike in the case of stable particles 
considered in Sec. |4.3| , the latter does not have a Gaussian form. Note the different time 
dependence of the different terms in the exponent in the integrand of the first integral in (|6^ : 
the terms proportional to a and b are multiplied hy t—tp, whereas the term oc F is multiplied 
by t. This is because the former two terms reflect the fact that the peak of the wave packet 
of Pi is located at the neutrino production point x = xp at the time t = tp, whereas the 
wave function of this particle exhibits an overall exponential suppression starting from its 
creation time t = 0. 

Calculating the integral in (|62|) , we find 



^ , , NpM.p(Q,K) r (p-P)^ r r b\[ / ib ^ \ 1 

(64) 

where erf(x) is the error function. The limiting cases of interest can now be obtained from 
the relevant expansions of this function, but it is actually easier to study them starting 
directly with the expression for Ji in eq. (|62D. 

Indeed, the integrand of Ii contains the exponential and oscillating factors, and therefore 
the integration domain that gives significant contribution to Ji is 

Itl < minj^, i| , (65) 

provided that the right hand side of this condition does not exceed tp] if it does, for negative 
t the domain is limited by |t| > tp. Let us note that, while the parameter b may vanish, 
b cannot, as it has a non-zero imaginary part r/2. Consider first the limit ^/a ^ \b\, 
y/atp ^ 1, which implies 

aeP>r/2, aeptp>l. (66) 



In this case one can set 6 ~ 6 and also extend the lower integration limit in eq. (|6^ ) to — oo, 
which gives Ji ^ a/tt/o exp[— 6^/4a]. Substituting this into eq. (|6ll) yields, up to the extra 
factor e"*^^/^^*^, the old result of eqs. ( ^^ and (^U]). If instead of the second condition in (|H^ ) 
one considers the opposite limit aeptp <^ 1, the lower integration limit in the last integral 
in (|62D can be set equal to zero. Since the error function goes to zero for small arguments, 
if follows that the result in this case is just 1/2 of that in the case a^ptp ^ 1. Thus, we 
conclude that in the limit a^p ^ r/2 the approximate conservation of mean energies is 
given by the same (in this case Gaussian) law as in the case of the stable neutrino source, 
with the same energy uncertainty cXep. This is an expected result. 
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Consider now the limit a/o ^ |&|, y/atp <^ 1, or 

cTeP <r/2, aeptp<l. (67) 

In this case one can neglect the term —at"^ in the exponent in the last integral in (|62D , which 
yields 

h i5,(p)-i?p-vp(p-P) + zr/2- ^ ^ 

This is the usual Lorentzian energy distribution factor corresponding to the decay of an 
unstable parent state. 

Thus, we conclude that in the case when r/2 ^ a^p the factor in $jp(i?j(p), p) that 
is responsible for the approximate conservation of mean energies in the production process 
is essentially the same as in the case of a stable neutrino source (Gaussian in the case we 
considered), whereas in the opposite limit, r/2 3> cXep, it is given by the Lorentzian energy 
distribution corresponding to the natural linewidth of the source r/2. In the intermedi- 
ate case r ~ (Tep, the energy- dependent factor in $jp(_Ej(p), p) is neither Gaussian nor 
Lorentzian, with the effective energy uncertainty being of the same order as OeP and F. 



5 Oscillation probabilities and the normalization con- 
ditions 



5.1 Normalization of oscillation probabilities in the QM approach 



In the QM wave packet approach to neutrino oscillations, one has to normalize the oscillation 
probability Pa/siL) by hand by requiring that it satisfy the unitarity constraint 

J2Pa^iL) = l. (69) 

This is an ad hoc procedure that is not properly justified within the QM formalism; how- 
ever, it seems to be unavoidable in that approach. In particular, one can readily make sure 
that the standard normalization of the neutrino wave packets = 1, or, equivalently, 

/ l/j(p)P'^^P/(27r)^ = 1, does not lead to the correct normalization of the oscillation prob- 
ability (that this is indeed the case can be easily verified by using Gaussian wave packets 
as an example). Moreover, as we will show below, no independent normalization of the 
produced and detected neutrino states can lead to the correct normalization of the oscillation 
probability in the QM wave packet approach. 

^°Note that the second condition in ( |67| ) cannot be relaxed: indeed, if one had aeP tp > 1, then r/2 ^ aeP 
would imply Ttp ^ 1, and the wave function of the parent particle would be exponentially suppressed by 
the time the neutrino is produced. 
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Let us now consider the unitarity constraints on the oscillation probabilities and their 
connection to the normalization conditions in more detail. As it is presented in eq. the 
unitarity condition simply reflects the fact that during the propagation between the source 
and detector neutrinos are neither destroyed nor (re) created: once a neutrino is produced, 
it can only change its flavour. Therefore, for a flxed initial flavour a, the probabilities 
Pa/3{L) of neutrino conversion to all flnal flavours (3 sum to unity. Note that the quantities 
Pap{L) by construction depend only on the distance L between the source and the detector 
and not on the propagation time T, which is integrated over (see eq. (|16D ) . Thus, the 
meaning of the unitarity condition (^) is that, once a neutrino is produced, the probability 
that it will be found at a given distance L from the source at some time between zero and 
inflnity is equal to one, provided that all flavour states are accounted for. Similarly, if one 
introduces Pq^(T) = J d^L P^piT, L), it will also satisfy a unitarity condition analogous to 
(p^). In other words, once the neutrino is created, the probability to flnd it (in any flavour 
state) at a flxed time T after its production somewhere in space is equal to one. 

The unitarity conditions can only be satisfled if the oscillation probabilities are properly 
normalized. It is important to note, however, that in a consistent formalism unitarity must 
be satisfled automatically rather than being imposed by hand. 

How about the un- integrated probability Poi^(T, L), should it satisfy a unitarity con- 
straint similar to (|69D? Obviously, for arbitrary L and T the answer to this question is 
negative. Indeed, unless |L — vT| < ax where v is the average group velocity of the neutrino 
wave packet and is its spatial length, the probabilities Pa,0[T, L) are vanishingly small for 
all a and /3. Thus, unitarity cannot be used to normalize the un- integrated probability. The 
normalization of Po^(T, L) can, however, be flxed differently: in the limit T — )■ 0, L — )■ 0, 
i.e. when the produced neutrino did not have time to evolve yet, the probability must satisfy 
the initial condition Pq/3(0,0) = Sap- Note that the above limit should be understood in 
the sense that L and T are small compared to the oscillation length but large compared to, 
respectively, the sizes of the spatial localization regions and time scales of the neutrino emis- 
sion and absorption processes. From eqs. (|T^ and (|T^) it then immediately follows that the 
amplitudes Aj(T, L) corresponding to neutrino mass eigenstates must satisfy Aj{0, 0) = e**^ 
where the real phase (p is the same for all j. By a rephasing of the momentum distribution 
functions of either emitted or detected neutrinos, this common phase can be eliminated, 
and one flnally gets ^^(0, 0) = 1. 

Now, let us check if this condition is fulflUed in the particular case of Gaussian wave 
packets. For the momentum distribution functions of the produced and detected neutrino 
states normalized according to eq. (^, i.e. having the form ([T0|), a straightforward calculation 

^^We neglect tiny probabilities of neutrino decay and absorption. 

^^For simplicity we assume here the neutrino emission to be spherically symmetric. Otherwise, Pq^(L) 
has to be defined as Pap{L) = J Pap{'L)dil/ATT. 
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yields 



For cTpp 7^ o"p£) and P 7^ P' both factors on the right hand side of the last equality are 
smaller than one, and therefore the condition .4.^(0, 0) = 1 is clearly violated. Moreover, the 
dependence of the result on the parameters of the produced and detected neutrino states 
does not factorize; this means that no independent normalization of these states can lead 
to the correct normalization of the amplitude, as pointed out above. 

It is actually quite easy to understand why this happens. The integral in eq. ([70|) 
is nothing but the overlap integral of the wave functions of the produced and detected 
neutrinos. If these wave functions are normalized to unity and the momentum (or energy) 
spectra of the emitted and detected states do not coincide, this overlap integral is always less 
than one. In reality, the spectra of the emitted and absorbed neutrino states are determined 
by the physical nature and experimental conditions of the neutrino production and detection 
processes, which are always different. This, in particular, means that a fraction of the 
produced neutrinos may simply not be detectable. For instance, if the threshold in the 
detection process (either the physical threshold or the one imposed by energy cuts of the 
detected events) is higher than the maximum energy of the emitted neutrino, no detection 
will be possible at all. Mathematically, the fact that the overlap integral (pO|) is always 
less than one is a consequence of the Schwarz inequality |(/,5')P < {f,f){g,g), where the 
equahty is only reached if / = const ■ g. 

Even if one adopts the unrealistic assumption /jp(p) = /jd(p) (which for Gaussian 
wave packets would mean P = P' and app = (TpD)-, this will not solve all the normalization 
problems of the QM wave packet approach. The condition ^^(0, 0) = 1 will be satisfied 
in this case; however, the physically observable oscillation probability Pq,^(L) defined in 
eq. (16) will still not be properly normalized, and the unitarity condition (|69|) will not be 



satisfied. Indeed, from eq. ( P^ ) it follows that unitarity requires that |T3 



j rfr|^,(T,L)|2 = l (71) 



for all j. Obviously, the fulfilment of the condition ^^(0, 0) = 1 does not enforce (fflD ; 
therefore in the QM formalism condition ( |7T1) has to be imposed by hand. 

It is not difficult to understand why yet another normalization problem arises: integra- 
tion over the time T should actually be considered as time averaging in the QM approach. 



and the integral on the right hand side of eq. (16) should be normalized by dividing it by 
the characteristic time AT which depends on time scales of both the neutrino production 
and detection processes. This follows from the fact that the amplitude Aj{T, L) is substan- 
tially different from zero only when |L — vT| < cr^., where the effective length of the wave 
packet (Jx is determined by both the neutrino production and detection processes, which 
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gives AT ~ ax/v. It is difficult to calculate the quantity AT precisely, and the simplest way 
out is just to impose the unitarity condition by hand, which yields the correct normalization 

Thus, in the QM approach there are two sources of the normalization problems: lack 
of the overlap of the wave functions of the produced and detected neutrino states and the 
necessity of integration over the neutrino propagation time T. We will show now how both 
these problems are naturally solved in the QFT-based formalism. 



5.2 Normalized oscillation probabilities in the QFT framework 
5.2.1 Generalities 

Let us start with recalling the operational definition of the neutrino oscillation probability. 
In a detection process that is sensitive to neutrinos of flavour /3, the detection rate is 

Vf = j dEj^{E)a^{E), (72) 

where o-p{E) is the detection cross section and j 13(E) is the energy density (spectrum) of 
the i>p flux at the detector. If a source at a distance L from the detector emits neutrinos of 
flavour a with the energy spectrum dr^°'^{E)/dE, the energy density of the z/^ flux at the 
detector is 

where PapiL, E) is the Va Vb oscillation probability for neutrinos of energy E, and we once 
again assumed for simplicity neutrino emission to be spherically symmetric. Substituting 
this into eq. ( [72D yields the rate of the overall neutrino production-propagation-detection 
process: 



The oscillation probability can now be extracted from the integrand of eq. ([7^) by dividing 
it by the neutrino emission spectrum, detection cross section and the geometrical factor 
l/47rL2: 

dT':l{E)/dE 

j^[dTr\E)/dE]a,{E) 

Note that an important ingredient of this argument is the assumption that at a fixed neutrino 
energy E the overall rate of the process factorizes into the production rate, propagation 



p^^iL, E) = . .,^p:d,J;,^. • (75) 



^•^Note that once the normalization condition ( [7l|) is enforced, one can demonstrate that the resuhing 
oscillation probability PapiL) is Lorentz invariant ||15|. This is not trivial in the QM approach because the 
QM formalism is not manifestly Lorentz covariant. 

^"'We omit the obvious factors of detection efficiency and energy resolution that are not relevant to our 
argument. 
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(oscillation) probability and detection cross section. Should such a factorization turn out 
to be impossible, the very notion of the oscillation probability would lose its sense, and one 
would have to deal instead with the overall rate of neutrino production, propagation and 
detection. 

Now let us come back to the QFT-based treatment of neutrino oscillations and try to 
cast the rate of the overall process in the form of eq. (0). To this end, we return to eq. (|2^ ) 
for the amplitude of the process but, unlike in the previous sections, perform in it the 
integration over 3-momentum before integrating over the energy variable p'^ = E. In doing 
so, we will make use of a result obtained by Grimus and Stockinger ||20|, which states that, 
for a large baseline L, positive A and a sufficiently smooth function ^/'(p). 



(76) 



whereas for A < the integral behaves as L^^. This result was obtained in in the limit 
L — )■ oo, but a careful examination of the derivation shows that its applicability condition is 
actually L ^ Pj/cTp, which we will assume to be satisfied. Applying ([76|) to eq. ( ^4]) yields 

tA^piTM = JdE^,piE,p,l)^,niE,p,\) 2E e-'^^^^^^^^ (77) 



where 



p, = ^E'^-m], \=^. (7J 



Next, we note that, just as the quantities $jp(i?j,p) depend on the index j only through 
the neutrino energy Ej^ the functions $jp(-E,Pjl) depend on the index j only through the 
neutrino momentum Py Therefore, to simplify the notation we will denote 

^jp{E,pj\) = <^p{E,p,\) , <^ME,Pji) = ^d{E,PjI) . (79) 

The overall probability of the neutrino production-propagation-detection process P^^(T, L) 
is the squared modulus of the amplitude ([T?]): 

P^°;(T,L) ^ |A^(T,L)|2 = J2u:,U^,U^,U;,A,iT,L)Al{T,L) . (80) 

We will actually need the integral of this probability over the time T (the reasons for 
integration over T will be discussed in Sees. |5.2.2| and 

X JdE^p{E,pjl)^n{E,Pj\)^*p{E,Pkm*D{E,Pkl) (2^)^ 6^^"^"^)^ . (81) 
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We use the tilde here to stress that P*^*(L) is not a probabihty but rather a time-integrated 
probabihty, which has the dimension of time. Note that (^Tj) contains an incoherent sum 
(integral) over contributions from different energy eigenstates. This means that only the 
amplitudes corresponding to the same neutrino energy interfere, as it is in the stationary 
case. This is related to the integration over time T and is a reflection of the fact that the 
time-integrated nonstationary probability is equivalent to the energy-integrated stationary 
probability |15|. It should also be stressed that, although the integration over E in eq. (|8T|) 
is formally performed over the interval (—00, 00) (recall that E coincides with the variable 
of eq. (p^)), the contribution of the unphysical region of negative energies is actually 
negligible and can be discarded. This is a consequences of the fact that $p,z) are sharply 
peaked at positive values of energy Ep^o, with the peak widths satisfying aep,eD ^ -E'p,d- 

To simplify the following consideration, we will once again assume that the neutrino 
production and detection processes are isotropic. In our approach this means that we have 
to average the quantities ^p{E,pjl) and ^DiE,pjl) over the directions of the incoming 
particles Pi and Di, which amounts to averaging over the directions of 1. We can therefore 
denote <^p^p,{E,pj) = J ^$p£)(i?,pjl) and drop 1 from the arguments of $p,d in eq. (|8T|) 
and all the subsequent expressions. Relaxing the isotropy assumption would complicate the 
analysis but would not change the final result for the probability of neutrino oscillations. 

The next step is to calculate the neutrino production and detection probabilities. As can 
be seen from eq. (|5), neutrinos in the intermediate state are considered in our framework 
as plane waves weighted with the factors $p,_d- We therefore describe the external particles 
by wave packets and neutrinos by plane waves. Apphcation of the standard rules of QFT 
then yields 

= E / "'^ l*f (^S-PjOr^SPj . (82) 

j 

The spectral density of emitted neutrino flux, dP^™'^(E)/dE, is obtained by removing the 
integration over energy on the right hand side of the last equality in (0). 

For the detection probability we obtain 

P''\E) = J2\U,k\'\^D{E,p,)\'^, (83) 

k 

where the normalization volume V comes from the plane-wave description of the incoming 
neutrino. Note that the expression for the production probability pp™"^ in eq. ( p2D does not 
contain the factor 1/V even though it is also calculated for plane-wave neutrinos. This is 

^^Note that we do not include the factor [2Ej{p)]^^ in tlic integration measure because it is aheady 
inchidcd in the definition of |<I>jp(i?, p)p, see eqs. ( p3| ) and (|25|). 
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because neutrinos are in the final state at production, and the calculation of their phase 
space volume involves integration over Vd^p. 



5.2.2 The case of continuous fluxes of incoming particles 

A direct inspection of the expressions for the probabilities of neutrino production and de- 
tection as well as of the probability of the overall production-propagation-detection process 
obtained above shows that they are independent of the total running time of the experiment 
t. This is because they were calculated for individual processes with single external wave 
packets of each type, and the microscopic production and detection time intervals were as- 
sumed to be centered at fixed instants of time tp and t^, respectively. On the other hand, in 
practice one is usually interested in the total probabilities for the processes to occur within 
a macroscopic time interval of length t or in interaction rates. Normally, the probabilities 
are proportional to t, while the rates are t-independent. In the wave packet approach, this 
can be achieved if we take into account that in realistic situations one often has to deal 
with continuous fluxes of incoming particles. (We will comment on the opposite case of 
stationary initial states in the next subsection.) 

Let us start with calculating the production rate and detection cross section in the case 
of steady fluxes of the incoming particles. Consider some interval of time Tq that is large 
compared to the time scales of the neutrino production and detection processes. Let the 
number of the projectile particles Pi entering the production region (e.g. a spherical region 
of radius axp around the point x = xp) during this interval be Np. Then the number of 
particles Pi entering the production region during the interval dtp is dNp^ = Np{dtp/TQ). 

If the production probability in the case of the individual process with single external 
wave packets is P^™'^, the probability of neutrino emission during the finite interval of time 
t (0 < t < To) is 16 

Vr'^it) = Np [ ^ Pf°^ = Np Pr'^^ . (84) 



„ u -1-0 

The integration is trivial because pp"^"'^ is actually independent of tp due to invariance with 
respect to time translations. Note that the probability V^^°'^ is proportional to t. We can 
therefore define the production rate in the usual way: 

Let us now consider the detection cross section. Normally, a cross section is defined for 
a single, fixed target particle. However, the initial-state particles Di in our treatment of the 



16 



If the flux of the incoming particles is not steady, the number of particles entering the production 



region over the time t is given by /q pp{tp)dtp, where pp{tp) is the distribution of these particles with 



respect to tp, normalized according to pp(tp)dtp = Np. The right hand side of the first equality in 

eq. (U) then has to be replaced by pp{tp)P^™^dtp. For a steady flux the distribution is uniform, i.e. 

pp{tp) = Np/Tq = const. 
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detection process are described by moving wave packets, which enter the detection region 
that is centered at a fixed point x = x^). Therefore, our treatment of neutrino detection 
should be similar to that of the production process. For the individual detection process 
with single external wave packets of each type the detection probability is given by eq. (pSf). 
Let us now assume that the number of the particles Di entering the detection region during 
the interval of time Tq is Nj^. Then we obtain for the time- dependent detection probability 
in the case of steady incoming fiuxes of Di and neutrinos 



Vf\t) = N, 



D 



To 



t 



and for the detection rate 



pdet 
^0 



dvrit) 



D 



pdet 

\3 



(86) 



7) 



dt Tq 

To obtain the detection cross section we have to divide this rate (more precisely, the sum- 
mand of the sum over k that enters into {^7\)) by the fiux of incoming neutrinos j^k = ^vkVuk 
where Ui^k is the number density of the detected Vk and v^k is their velocity. With our nor- 
malization (one particle in the normalization volume) we have Uyk = 1 1^1 and from eqs. (|8 
and (|87D we finally obtain 



<jp{E) 



\^o{E,pk)?-. 

Pk 



Now we proceed to the calculation of the rate of the overall production-propagation- 
detection process. Since we want to calculate this quantity not for a single process with 
individual wave packets of external particles but for steady fiuxes of incoming Pi and Di, 
we have to integrate the T-dependent probability of the single process P^^(T,L) given 
by eq. (|80|) over both tp and tn. Proceeding in the same way as before, we find for the 
probability of the process for steady fiuxes of the incoming particles 



NpNd 
-'o 



dt 



D 



dtpnt{T,L) 



^9) 



Introducing the new integration variables T = {tp + tp,)/2 and T = tp) — tp, we obtain 



-'o '-"'0 



dTPl^'{T,L){t-T) 



dTP'^f{T,L){t + T) 



NpNo 
NpNo 

-^0 



t j dTP'°^{T,L)- / dTTP'J^;{T,L) 
th{t) - hit) + hit) 



dTTP'^'iT,L) 



(90) 



It can be readily shown that in the limit of large t (much larger than the time scales of 
the neutrino production and detection processes) the integral h coincides with the quantity 
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P^p{L) defined in eq. (^TJ), whereas I2 and /s give negligible contributions (see Appendix B). 
Therefore for large t eq. (0) can be rewritten as 

V:;{t,L) = ^^tPl^^{L). (91) 
The rate of the overall process is then 

C(i) = ^^^ = A'pA'«^. (92) 



5.2.3 The case of stationary initial states 

If the initial state particles Pi and/or Di are in stationary states rather than being described 
by moving wave packets, the above consideration has to be slightly modified. Consider, for 
example, neutrino production in decays of unstable particles Pi bound in a solid. Let pp{tp) 
be the probability distribution function for the decay times of the parent particles in the 
source, and let Np be defined by the condition pp{tp)dtp = Np. If Tq is short compared 
to the lifetime F"^ of Pi, one has pp{tp) = const = Np/Tq. If, on the contrary, Tq > F"^, 
pp{tp) will usually have an exponential form.^^ The neutrino production probability is then 
again given by eq. (|84D, just as in the case of a continuous flux of incoming particles Pi. 

The situation with bound-state stationary particles Di in the detector can be considered 
quite similarly. One can assume Di to be stable. If the source creates a steady flux of 
neutrinos, then for the ensemble of Di in the detector the distribution puitp,) oi the detection 
times tp) is uniform and is given by Nd/Tq, where Np, is defined by the normalization 
condition Jq" p 0(^0) dtp, = Np. The detection probability is then given by eq. (^). 

Thus, with these re- interpretations of Np and A'^^, the expressions for the neutrino 
production and detection probabilities and rates and the detection cross section obtained 
in the previous subsection remain valid in the case of stationary initial states as well. 

We are now in a position to obtain the normalized oscillation probability. 



5.2.4 The oscillation probability in the QFT approach 

In the case when the rate of the overall production-propagation-detection process at a fixed 



neutrino energy factorizes, the oscillation probability should be obtainable from eq. (75). 



^^For the reader willing to check the dimensions of our expressions, we note that from the definitions of 
<f>f,_D it follows that they have dimension m""^/^. It is then easy to see that the amplitude ( [77|) and the 
probabilities (|8ll)-(^4|), (|86|), ( p9| ) and ( [qTI) are dimensionless, the cross section ( ^8|) has the dimension of 
squared length (or m^^), and the rates (|85|), ( p7| ) and ( |92| ) have the dimension of inverse time (or m), as 
they should. 

^^An exception is the case where Pi in the source are continuously replenished. In this case, the function 
Ppitp) will depend on the time dynamics of the production of these particles. 
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Substituting eqs. (H) and (||) (with P^f and PP'°'^ defined in eqs. (|T|) and (|D) and 



eq. into this relation, we find 

(93) 

The quotation marks here are to remind us that we yet have to prove that this quantity can 
indeed be interpreted as the oscillation probability. 

The alert reader has probably noticed that, while the integral in ( [7^ ) is taken over the 
energies of different neutrinos in the neutrino flux, the integration in (^) is performed over 
the energy distribution within the wave packet of an individual neutrino. This, however, does 
not invalidate our argument leading to eq. (^3|) . The reason for this is that the following two 
situations are known to be experimentally indistinguishable (a) a flux of neutrinos 

described by identical wave packets, each with an energy spread f{E), and (b) a flux of 
neutrinos, each with a sharp energy, with the overall energy distribution (f){E) = \ f{E)\'^. 

Let us now examine expression (^3]). First, we note that when the produced neutrinos 
are either ultra-relativistic or quasi-degenerate in mass, i.e. when \pj—pk\ Pj,Pk, the prob- 
abilities of emission of different neutrino mass eigenstates are characterized by essentially 
the same transition matrix elements and the same phase space volumes, i.e. to a very good 
accuracy these probabilities only differ by the factors Therefore, one can replace the 

factors ^p{E,pj) in eq. (^) by the value ^p{E,p) calculated at the average momentum 
p, and also replace the factors pj by p, and pull them out of the sum. A similar argument 



applies to the detection cross section. As a result, in the denominator of ( p3D we can replace 



J2\Uaj\'\<^piE,p,)\'p,^\<^p{E,p)\'pY,\Ua,\'=\'^p{E,p)\'p, 

j j 

If/ftf \<l>D{E,p,)\'p-,' -> \<l>n{E,p)\'p-' \Up,\' = \<l>n{E,p)\'p-' , (94) 

k k 

where in the last equalities we have used unitarity of the leptonic mixing matrix. Note that 
such a procedure cannot in general be applied to the numerator of eq. (P^. Indeed, the latter 
contains interference terms proportional to products of the functions $p taken at different 
momenta, pj and pk, where pj was deflned in ([781). The momentum distribution functions $p 
are all peaked at the same momentum P, therefore, if \pj —pk\ exceeds the width of the peak 
CTpp, the interference terms will be strongly suppressed. The same argument applies to the 
detected neutrino states $pi(£',pj) and $£>(£', p^), which also enter into these interference 
terms and which are peaked at momentum P' and have widths apD- For ultra-relativistic 
or quasi-degenerate neutrinos (i.e. when \pj — pk\ <^ p), one has \pj — Pk\ — Am^^/2p; if 



^^It should be stressed that the mean momentum p is defined here as an average over different mass 
eigenstates of the momenta pj = {E"^ — m|)^/^ taken at the same fixed value of energy E. It is therefore 
different from the mean momentum P of the individual wave packets, introduced earlier, for which the 
average was taken over the spread of momenta (or energies) within the wave packet. 
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Am^^/2p > (jp where is the effective momentum uncertainty dominated by the smallest 
between (jpp and apo-, the interference terms in (and therefore neutrino oscillations) will 
be strongly suppressed. Physically, this can be traced to the lack of coherence at neutrino 
production and/or detection. It can be shown that production or detection decoherence is 
equivalent to the lack of localization of, respectively, the production or detection process 
|jl|, |l5l. If, on the contrary, 

^ « , (95) 

i.e. the difference \pj—pk \ is small compared to the effective width of the neutrino momentum 
distribution, the production/detection coherence condition is satisfied, and the oscillations 
can be observed. Since the effective momentum uncertainty ap is actually dominated by 
the energy uncertainty (Tg (see Sec. ^73|) , condition (^) is equivalent to the one in eq. (|60|). 

By comparing eqs. (|8ll) and (|93|) , it is easy to see that in the approximation (Q), 
when the spectral density of the production probability and the detection cross section are 
independent of the elements of the leptonic mixing matrix, they can be factored out of the 
sum over j and k in eq. ([HT|), which leads to a factorization of the form ( [7^ ) . This means 
that the oscillation probability Pai3{E, L) can be defined as a sensible quantity and is given 
by 

"^^^^^'""^^ \ME,pW\^niE,pW • ^''^ 

With the help of eq. (0), it is easy to make sure that this expression automatically satisfies 
the unitarity condition (^), i.e. is properly normalized. Thus, the QFT-based approach 
allows one to identify the conditions under which Pq,^(L, E) can be sensibly defined, and also 
gives the correctly normalized expression for this probability. The condition for the existence 
of well-defined oscillation probabilities is that neutrinos are either ultra-relativistic or nearly 
degenerate in mass. Note that the factorization of the integrand of (0) according to (|7^) 
does not by itself mean that the oscillation probability ( p6|) is production and detection 



independent; for this, one would also have to require that condition (pSQ be fulfilled 
If it is satisfied, all the momenta Pj,Pk in the interference terms are sufficiently close to 
each other, and one can replace ^p{E,pj) — > ^p{E,p), ^pi{E,pj) ^pi{E,p) also in the 
numerators of eqs. (|93D and (p6|). As a result, the factors |<l'p(-E',p)p |$£i(-E',pfc)p in the 
numerators and denominators of these equations cancel, and one obtains 

P.,(L, E) = J2,/*»^Upj^o.kU;, e^(^^-^'=)^ . (97) 



While condition (|95| ) ensures the production/detection coherence (locaUzation), it says nothing about 
another possible source of decoherence - separation of neutrino wave packets at long enough distances 
L > Lcoh due to the difference of the group velocities of different neutrino mass eigenstates. This is related 
to the fact that a fixed neutrino energy corresponds to the stationary situation, when the coherence length 
^coh oo. The finite coherence length is recovered upon the integration over energy in eq. (pw B. 
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Since for ultra-relativistic or quasi-degenerate neutrinos pj — pk — —^m?-^,/2p, this is just 
the standard formula for the probabihty of neutrino oscillations in vacuum. 

The above QFT-based considerations also allow one to shed some light on the meaning 
of the normalization condition imposed on the oscillation probability in the QM wave packet 
approach, which looks rather arbitrary within the QM framework. As was pointed out in 
Sec. the QM and QFT approaches can be matched if the QM quantities fjp and fjo 
are identified with the QFT functions <l>jp(i?j,p) and $*^(_Ej,p), respectively. The latter, 
however, bear information not only on the properties of the emitted and absorbed neutrinos, 
but also on the production and detection processes. The QM normalization procedure that 
is tailored to obtain an expression for the oscillation probability satisfying the unitarity 
condition can then be easily seen to be equivalent, in the limit (CT), to the division of the 



overall rate of the process by the production rate and detection cross section, as in eq. (96). 



6 Some additional comments 

In this section we will comment on some issues pertaining to the description of neutrino 
oscillations in the external and intermediate wave packet approaches that were not discussed 
or were only briefly mentioned above. We will also show how one can relax some assumptions 
usually adopted in the QM and QFT approaches. 

1. Unequal mean momenta of the produced and detected neutrino states 

In the vast majority of derivations of the neutrino oscillation probability within the 
QM wave packet framework, it was assumed that the mean momenta of the produced and 
detected neutrino states, P and P', coincide. There are, however, no reasons for this to 
be the case. Indeed, the mean momenta of the emitted and absorbed neutrino states are 
determined by the kinematics and experimental conditions of the neutrino production and 
detection processes, respectively, and those are in general different. Let us now examine the 
consequences of P 7^ P' using, as before, the case of Gaussian wave packets as an example. 

It was shown in Sec. |4.3| that the neutrino wave packets derived in QFT can be cast into 
the form usually adopted for neutrino wave packets in the QM approach if one uses for the 
latter the effective (in general anisotropic) momentum uncertainties, shifted mean momenta 
and modified normalization factors. Since the expressions for the QM wave packets look 
simpler, we will study the implications of unequal mean momenta of the produced and 
detected neutrino states within the QM formalism. We will also ignore possible anisotropy 
of the neutrino momentum uncertainties. Taking it into account would just complicate the 
calculations without changing the essence of the result. 

We will now calculate the amplitude ([T5|) corresponding to the emission and absorp- 
tion of the neutrino mass eigenstate Uj. To do so, we expand the energy -E'j(p) around a 

^^The only exceptions we are aware of are refs. |27|] . 
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momentum po which we do not specify for now: 



£:,(p)~£:,(po)+v,(p-po). (9^ 



An important point is that, strictly speaking, using such an expansion in the integral in ([T5|) 
is only justified if P and P' are not too far from each other and if po = 0(P, P'), which we 
will assume. As we shall see, under these conditions the final result will be insensitive to 
the choice of the expansion point po- In eq. ( pS]) we have neglected higher order terms in 
the expansion; this corresponds to neglecting the spreading of the neutrino wave packets, 
which is unimportant for our argument. 

The Gaussian wave packets for the produced neutrino states are given by eq. (PUj), and 
the detected wave packets have a similar form, with app and P replaced by and P', 
respectively. Substituting these expressions and expansion ( P5| ) into eq. (^) yields, after a 
simple Gaussian integration. 



^^.(T, L) ^ I ^fp^^ \ e *(ip+^^z,)e-^.(P)^+'Pi^ e"^^ . (99) 



3/2 



CTpP + (Jpj^ 



Here 



__Pal^ + P^ 1 _ 1 1 _ 1 

(Jpp -h (Tp^ (Tp (Tpp (7p£, Za„ 



and we have used the relation -E'j(po) + Vj(p — po) — Ej{p) which follows from (^. As can 
be seen from eq. (0), the dependence on the expansion point po has completely disappeared 
from the amplitude. 

The meaning of the obtained result is quite transparent: as usual, the amplitude 
Aj{T,'L) contains a normalization factor, a plane wave factor calculated at some mean 
momentum (in this case p) and the envelope factor exp[— (L — VjT)^/4(T^]. However, on top 
of this, it contains the extra factor exp[— (P — P')^/4(crpp + (jp£,)], which is a reflection of the 
approximate conservation of the mean neutrino momentum. It suppresses the amplitude 
unless the difference of the mean momenta of the produced and detected neutrino states is 
small compared to the effective total momentum width (o'pp + CTpp,)^/^. Note that this width 
is dominated by the larger of app and <TpD, unlike the momentum width ap defined in (|100|) , 
which is dominated by the smaller of them. 

2. Non-central collisions 

In our discussion of the external wave packet formalism, as well as in all other treatments 
of this topic we are aware of, it was assumed that the peaks of the wave packets of all the 
external particles participating in neutrino production (or detection) meet at the same 
space-time point. In other words, it was assumed that at production these peaks are all 
located at the point x = xp at the same time t = tp, whereas at detection the peaks of the 
wave packets of the external particles are all located at x = X£) at the same time t = tp,. Is 
this always true and how crucial is this assumption? 
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While in some situations, such as neutrino production in decays of unstable particles, 
such an assumption may indeed be justified, this is not in general so when neutrinos are 
born in scattering processes, where the collisions of the wave packets in the initial state 
may be non-central. Let us consider non-central collisions assuming, as before, only two 
external particles at production and two at neutrino detection (the generalization to an 
arbitrary number of external particles is straightforward). One obvious consequence of this 
is that the neutrino production and detection are only possible if the minimum distances 
between the peaks of the participating wave packets do not exceed significantly the sizes 
of, correspondingly, the production and detection regions. Consider the production process. 
We shall now assume that the wave packets of the external particles Pi and Pf are given 
by expressions of the type (^, with the phase factors in the integrand being ^-i-ii^-^'^) 
and e-«^(^-^i>)^ respectively. This means that the peak of the wave packet describing Pj is 
located at x = at the time t = ta, and the peak of the wave packet of Pf is located 
at x = Xf, at the time t = tf,. We will consider the positions of these peaks at the same 
time, i.e. we will take ta = tb- It is natural to choose tp to be the value of this common 
time corresponding to the minimum distance between the two peaks, i.e. ta = tb = tp, 
|xa — Xb| = min{|xa(t) — Xb(t)|}.^^ The coordinate xp can then be chosen to lie anywhere 
between Xq and x^ on the line connecting them; in particular, one can choose xp = Xq or 
Xp = Xf,. Neutrino detection can be considered quite similarly. 

One can now repeat the calculations presented in Sec. 0, arriving at the transition am- 
plitude that can again be written in the form (|2^) . However, the expressions for $jp(p°, p) 
and ^joip^, p) in eq. ( pS]) have to be modified: their integrands should be multiplied, respec- 



tively, by e«q(xp-xa)-«k(xp-x6) ^^^^^ gjq'(xp-xc)-ik'(xp~xd)^ where Xc and x^ are the positions of 
the peaks of the wave packets representing Di and -D/ at detection at the time to when the 
distance between these two peaks reaches its minimum. This can be taken into account by 
redefining the momentum distribution functions of the wave packets according to 

/p.(q,Q) ^ /«(q,Q)e^^(''--^"), /p/(k,K) -> fpfik,K)e^^--~->'\ 

/D.(q',Q') -> /D.(q^Q')e^'^'(''--''=^ fof {k',K') ^ /z,/(k', K')e^'''(''--^^) . (101) 

The newly defined momentum distribution functions share a crucial feature with the old 
ones: they decrease rapidly when the deviation of the corresponding momenta from their 
peak values exceed the relevant momentum uncertainties (i.e. the widths of the peaks) 
o'pPi, o'ppf, (JpDi or (TpDf- In addition, the new momentum distributions do not exhibit 
fast oscillations when the momentum variations are smaller than, or of the order of, the 
corresponding widths of the peaks. Indeed, consider neutrino production. Our choice of xp 
implies that |xp — Xa| < |xa — X(,| and |xp — x;,| < |xa — Xf,|. On the other hand, as we 
mentioned above, the production is only possible when the distance between the peaks of 
the wave packets of the external particles |xa — x;,| is smaller than, or of the order of, the size 



^^This is easily generalized to the case of more than two external particles. For instance, for three particles 
with peaks at Xq, x^ and Xc the time tp can be chosen to correspond to the minimum of (x^ — x;,)^ + (x^ — 

Xc)2 + (Xfc - Xc)2. 
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of the production region, which is of the order of [max{(TpPj, (Tpp/}]"^. A similar argument 
apphes to neutrino detection. Therefore the variation of the exponents of the exponential 
factors in eq. (|101|) is < 1 and these factors do not undergo fast oscillations across the peaks 



of the momentum distributions. 

From the above considerations it follows that the properties of the redefined momentum 
distributions are essentially the same as those of the old ones. All the results of the present 
paper therefore apply to the case of non-central collisions as well, if one substitutes the 
original momentum distribution functions by those redefined according to eq. ( |10ip . 



7 Discussion and summary 

In this paper we have compared the quantum mechanical approach to neutrino oscillations, 
where neutrinos are described by wave packets, with the quantum field theoretical method, 
where they are represented by propagators connecting the neutrino production and detection 
vertices in a Feynman diagram, whereas the external particles are described by wave packets 
in order to localize the process in space and time. We have shown how the neutrino wave 
packets underlying the QM approach can be derived in QFT by comparing the QM and QFT 
expressions for the transition amplitude. Equivalently, the wave packet representing the 
emitted neutrino can be obtained as the convolution of the neutrino source (the production 
amplitude) with the retarded neutrino propagator, in accord with the well known result of 
QFT. Quite analogously, the wave packet of the detected neutrino can be obtained as the 
convolution of the neutrino detection amplitude with the advanced neutrino propagator, 
with the result then taken at the time corresponding to neutrino detection. 

We have studied the general properties of QFT-derived wave packets representing the 
produced neutrino states and demonstrated that the wave packets of mass eigenstates Vj 
depend on the index j only through the neutrino energy Ej, and that in the momentum 
representation they are given by the production amplitude multiplied by "smeared delta 
functions" describing approximate conservation of mean energy and mean momentum at 
production. The widths of these "smeared delta functions" are determined by the largest 
among the corresponding widths of the external particles involved in the neutrino production 
process. Similar conclusions apply to the wave packets of the detected neutrino states. 

We also identified the conditions under which general neutrino wave packets can be 
approximated by Gaussian ones. Using Gaussian wave packets as an example, we then 
demonstrated that the neutrino wave packets derived in QFT can be cast into the form 
they are usually assumed to have in the QM formalism, provided that 

(i) The momentum uncertainty of the QM approach is replaced by the effective one, 
which depends not only on the true momentum uncertainty at production (or detection), 
but also on the corresponding energy uncertainty, as well as on the neutrino velocity and the 
effective velocity of the neutrino production (or detection) region. Moreover, these momen- 
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turn uncertainties are different in different directions, i.e. are anisotropic. The longitudinal 
effective momentum uncertainties of the produced and detected neutrino states are domi- 
nated by the energy uncertainties characterizing, respectively, the neutrino production and 
detection processes, whereas the transverse effective momentum uncertainties coincide with 
the corresponding true momentum uncertainties. 

(ii) The mean momentum of the neutrino state is shifted from its naively expected value. 

(iii) The wave packets of different mass eigenstates acquire (in general different) extra 
overall factors. 

Thus, the simplistic QM wave packet approach may need QFT-motivated modifications; 
however, once these modifications have been done, one can still work within the QM frame- 
work without losing any essential physical content. 

We have also studied the energy uncertainties characterizing the neutrino wave packets 
in the case of unstable neutrino sources and have shown that in general these uncertainties 
depend both on the decay rate of the parent particle and on the inverse time scale of the 
overlap of the wave packets of the external particles. The neutrino energy uncertainty is 
dominated by the larger of the two. 

In the last part of the paper, we have discussed in detail the normalization of the 
QM and QFT expressions for the oscillation probability Pap{L). We have seen that in 
the QM framework PajS^L) has to be normalized by hand in order to fulfill the unitarity 
relation ^^-P«^(i^) = 1. There are two reasons why this ad hoc procedure is unavoidable 
in the QM approach. First, as we have demonstrated, no independent normalization of the 
produced and detected neutrino states can lead to the correct normalization of the oscillation 
probability because the overlap integral of these states is always smaller than unity in the 
realistic case when their momentum distributions are different. The overlap integral depends 
on the characteristics of the produced and detected states in a non-factorizable way, and so 
the problem cannot be cured by just modifying the normalization of these states. Second, 
the QM formalism involves an integration over the unobserved difference of the neutrino 
detection and production times T = to — tp, which leads to yet another undefined factor 
- the time interval by which one has to divide the result in order to recover the correct 
dimension of the oscillation probability. In the QM method both these problems are solved 
by imposing unitarity of the oscillation probability by hand. 

We have demonstrated how the QFT approach avoids all the normalization problems of 
the QM formalism and naturally leads to the correctly normalized oscillation probability that 
automatically satisfies the unitarity condition, provided that neutrinos are ultra-relativistic 
or quasi-degenerate in mass. If this requirement is not fulfilled, the interaction rate cannot 
be factorized into the production rate, propagation (oscillation) probability and detection 
cross section, so that the oscillation probability is undefined. In that case one would have to 
deal instead with the overall rate of the neutrino production-propagation-detection process. 

The QFT approach also allows one to understand the physical meaning of the QM nor- 
malization recipe: By imposing unitarity by hand one implicitly rids the calculated transi- 
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tion probability of the probabilities of neutrino production and detection, thus extracting 
the sought oscillation probability. 

A comment on the the integration over T is in order. Such an integration is involved in 
both the QM and QFT approaches to neutrino oscillations. In the QM framework, it has 
to be introduced to account for the fact that the neutrino's time of flight is not measured 
(or at least not measured accurately enough) in realistic experiments. At the same time, 
in our QFT treatment of neutrino oscillations, it emerges naturally from the observation 
that in real situations one has to deal with continuous fluxes of incoming particles (or with 
ensembles of neutrino emitters and absorbers in the case of bound stationary initial states) 
rather than with individual acts of neutrino production, propagation and detection, in which 
single wave packets of the external particles of each type are involved. Can one still sensibly 
define an unintegrated oscillation probability Pq,^(T, L) for such a single act? 

We will argue now that the probability Pq,^(T, L) is not a useful quantity since it is 
unmeasurable (or almost unmeasurable) . An important point here is that in practice both L 
and T can only be measured with some accuracy. If we consider the T-integrated probability 
Pap{L) which depends (in addition to neutrino energy) only on L, then this quantity is well 
defined only if the error AL in the measurement of L is small compared to the variations 
of L over which the probability changes significantly. This means that this error must be 
small compared to the neutrino oscillation length lose-, and this condition is normally easily 
satisfied. If one considers the case when the time T is measured, whereas the distance L is 
not (even though it is hard to imagine how such a situation could be realized in practice), 
then one would have to integrate P^^{T, L) over L, and the resulting probability would be a 
function of T. Then the situation would be similar - the error AT in the determination of T 
would have to be small compared to the change of T over which Pq,^(T) varies significantly, 
which is ~ (losc/v) with v the neutrino velocity. 

The situation would be quite different if one considers the unintegrated probability 
Pa/^iT, L) which depends on both T and L - in this case the requirements on AL and AT 
would be by far more demanding. Indeed, Pap{T,L) is substantially different from zero 
only when \L — vT\ < ax-, where ax is the spatial length of the neutrino wave packet. This 
means that for uncorrelated variations of L and T this probability varies significantly when L 
changes by ax (which can be extremely small) and T changes by (ax/v). The latter quantity 
is essentially given by the largest between the time scales of the neutrino production and 
detection processes. Thus, the un- integrated probability PapiT.L) can only be accurately 
measured if the distance L is measured with an accuracy better than the length of the 
neutrino wave packet and simultaneously the time between neutrino emission and detection 
is measured with an error that is small compared to the duration of these processes. Such 
a possibility appears rather unrealistic at the very least. 

It follows from our results that there are both intricate relations and important differ- 
ences between the QM and QFT approaches to neutrino oscillations. In the following table, 
we compare the main features of these two approaches. 
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QM approach 



QFT approach 



1. Simple and transparent description of neu- 
trino oscillations. Neutrino production and 
detection processes are not properly taken 
into account. Simplified description of neu- 
trino energy and momentum uncertainties. 

2. Produced and detected neutrino states are 
flavour eigenstates, defined according to 

3 

3. Mass eigenstates composing neutrino flavour 
eigenstates are described by wave packets 
whose form is postulated rather than derived 
and the parameters (momentum uncertain- 
ties) are estimated from the properties of the 
production and detection processes. Deriva- 
tion of neutrino wave functions is not pos- 
sible since they depend on the dynamics of 
neutrino production and detection, and par- 
ticle creation and annihilation cannot be de- 
scribed in QM. 

4. The oscillation amplitude is obtained by 
evolving the produced neutrino state in time 
and then projecting it onto the detected 
state. 



Most complete description of neutrino pro- 
duction, propagation and detection. Accu- 
rate treatment of neutrino energy and mo- 
mentum uncertainties. The formalism is 
more complicated than that of the QM ap- 
proach. 

Only mass eigenstates are considered. (In 
fact, defining flavour eigenstates in QFT 
poses great difficulties because they do not 
form a physically meaningful Fock space 




Because neutrinos are only in the intermedi- 
ate states, their wave functions are not nec- 
essary for the formalism (but can be derived 
from the production and detection ampli- 
tudes according to well defined rules). Wave 
functions of the external particles accom- 
panying neutrino production and detection 
have to be known. Undetected external par- 
ticles can be described by plane waves. 



The amplitude of the combined process of 
neutrino production, propagation, and de- 
tection is computed according to the Feyn- 
man rules. Time evolution of the neutrino 
states in QM corresponds to their on-shell 
propagation in QFT. Projection in QM cor- 
responds to the integration over the momen- 
tum of the intermediate neutrino in QFT. 

The oscillation probability Pap{L) that is 
properly normalized and satisfies the unitar- 
ity constraint is automatically obtained from 
the formalism in the case when neutrinos are 
ultra-relativistic or quasi-degenerate. Oth- 
erwise Pai3{L) is undefined. 



5. The oscillation probability has to be normal- 
ized by hand by imposing the unitarity con- 
dition. (The physical meaning and justifica- 
tion of this normalization procedure is eluci- 
dated in QFT.) 
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In summary, we have explicated the close relation between the quantum mechanical and 
quantum field theoretical approaches to neutrino oscillations and have shown how QFT - 
apart from providing expressions for oscillation probabilities and event rates in its own right 
- can be used to derive the input parameters required for the QM approach and to elucidate 
some QM procedures which were not properly justified or fully understood within that 
approach. We have also clarified several subtle points regarding neutrinos from unstable 
sources, the case of unequal mean momenta of the produced and detected neutrino states 
and the normalization of the oscillation probability. 
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Appendix A: Matching the QFT and QM wave packets 
in the 3-dimensional case 

We present here the full solutions of eqs. (0) and (|57D. 

The quantities a^'' and 6'^ depend on the components of two linearly independent vectors, 
Vj — vp and vj. We therefore seek the solution of eq. ( |56D in the form 



(Al) 



Substituting this into eq. (^G)) and comparing the coefficients of {vj —Vp) 
sides of the equality, one finds 



and of v\ on both 



Co = 



(E,-Ep)[Ap + ^(l-v|)] 



(A2) 



[Ap + - v|)][Ap + (v, - vp)2 + ^] + 



[{Wj - Vp)Vjf 



Ej - Ep [(vj - vp)vj] 



(A3) 



[Ap + ^(l-v|)]- 



The parameter 7^ can then be obtained from eq. (|57D, which gives 




(A4) 
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Substituting here the expressions for cq and ci, one finds 

„ {E, - Epf [Ap + ^(l-v|)][Ap + ^] 

4a,V [Ap + ^(1 - v|)][Ap + (v, - vp)2 + ^] + ^[(v, - vp)v,f 

(A5) 

The above expressions simphfy significantly if one makes use of the fact that 
is small (note that for relativistic neutrinos the above formulas contain one more small 
parameter, (1 — v|)). In the leading order in the effective momentum uncertainties 

that are related to the eigenvalues of the matrix a^' are given by eq. (|58| ) of Sec. ( [4.3| ), 
whereas the shift of the mean momentum and the parameters 7^ are given by eq. (|59D . 



Appendix B: Calculation of the integrals /i, 1^ and /g 



Here we calculate the integrals Ji, I2 and /a that enter into eq. ( pO]) in the limit of large t {t 
much larger than the time scales of the neutrino production and detection processes). For 
generality, we will not assume isotropy of neutrino emission and detection here. In the limit 
of large t the integral 1\ is given by 

h = ^^2Y.^iP,P.k^h r dT jdE l'dE'<^p{E,p,l)<l>oiE,p,l) 

x^*p{E',p'^l)^}j{E',p[\) {AEE')e-'^^-^'^^+'^P^-P'>^^^ , (Bl) 

where we use the shorthand notation p'^ = pk{E'). This quantity coincides with the time- 
integrated probability P^fiL) defined in eq. (|T]) since dT e''^'^-^")^ = 2'k5{E - E'). 

Let us now show that in the limit of large t the contributions of the integrals I2 and J3 
to eq. (pop are negligible. We first notice that is essentially zero because so is P^"^(T,L) 
for T < 0. From = it follows that we can replace I2 — I3 by I2 + I3, which turns out 
to be easier to calculate. For this quantity eq. ( ]90| ) yields /2 + /s = J2j k ^aj^pj^ak^pk ^jk 
with 

Bjk = ^Jr^/ dTT jdE jdE'<l>p{E,p,\)<l>n{E,p,\) 

x^*p{E',p',,r)^}j{E\p'i,l) (4EE')e-*(^-^')^+'(P^-P'^)^. (B2) 



Note that the integral here is similar to that in eq. (pTD, but, unlike the latter, contains 
J^^dTTe~'^^^~^^^ = —2Txi5'{E — E'), where the prime stands for differentiation with 
respect to E'. Because the integral over E' of f{E')6'{E — E') yields — f'{E), one obtains 



_ 2m 



r°° d 

j dEC,{E)—Cl{E), (B3) 
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where Cj{E) = ^p{E,pjl)^D{E,Pji) {2E) e-'^^+'P^^^^^. Integrating in (^) by parts we 
find 

= ^^mE)Cl{E)\-^ - IjECm — . (B4) 

The first term on the right hand side vanishes because so do the functions $p,z) at — t- ±oo. 
Therefore eq. ( p^) means Bjk = ~Blj. On the other hand, from the definition ( p^) of Bjk 
it follows that Bjk = B^-. Hence, Bjk = and I2 + 13 vanishes. 
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